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In this paper we experimentally and theoretically investigate laser cooling of Strontium 88 atoms 
in one dimensional optical molasses. In our case, since the optical cooling dipole transition involves 
a Jg = groundstate, no Sisyphus-type mechanisms can occur. We are thus able to test quantita- 
tively the predictions of the Doppler-cooling theory. We have found, in agreement with other similar 
experiments, that the measured temperatures are systematically larger than the theoretical predic- 
tions. We quantitatively interpret this discrepancy by taking into consideration the extra-heating 
mechanism induced by transverse spatial intensity fluctuations of the optical molasses. Experimen- 
tal data are in good agreement with Monte-Carlo simulations of our theoretical model. We thus 
confirm the important role played by intensity fluctuations in the dynamics of cooling and for the 
steady-state regime. 

PACS numbers; PACS; 32.80.Pj 

I. INTRODUCTION 

Initiated in the mid-seventies, laser cooling, trapping and manipulation of atoms rapidly became a very successful 
field of research Jj , culminating twenty years later with the observation of Bose-Einstein condensation of alkali atoms 
. This success story begun with the seminal Doppler cooling theory which was designed for two-level systems 0, 0] ■ 
However, the first reliable tests of the Doppler theory could only be performed at the end of the eighties. As it 
turned out, the measured temperatures in the experiments were well below the predicted Doppler values This 
surprising results indicated that another much more efficient cooling mechanism was at work in the experiments. This 
mechanism, now known as Sisyphus cooling, was identified soon after |6|. The key point was to understand that 
internal groundstate degeneracies were opening the way to a new cooling mechanism based on optical pumping in the 
presence of polarization gradients. This cooling mechanism is now the basic ingredient for most experiments in the 
field and thus gave the Doppler theory a more academic status. 

Recent laser technological advances however opened the way to cooling and trapping experiments with earth-alcaline 
and rare-earth atoms (Calcium Ca, magnesium Mg, strontium Sr and Ytterbium Yb). As these atoms exhibit a zero 
spin groundstate, Sisyphus cooling is absent. Interest in testing the Doppler theory to gain better control on achievable 
temperatures has thus been renewed. To our knowledge, all experiments on laser cooling of these atoms 0, 0, S IfOt 
always reported much larger temperatures than the Doppler theory prediction. This is a strong clue in favor of 
an extra-heating mechanism still to be clearly identified. Approaching the Doppler limit in magneto-optical traps 
(MOT) operating with earth-alcaline atoms seems a challenge and a better understanding of this heating mechanism 
is required for future experimental improvements . To our knowledge only one explanation has been published to 
understand these high temperatures The argumentation is based on heating induced by inelastic collisions. 

However, as the authors themselves show, this heating mechanism sets in only for atomic densities much higher than 
for those measured in typical MOT (10^ — lO^^atoms/cm'^). For this reason this explanation fails to understand the 
results presented in this paper and also in all previously cited publications. 

Heating mechanisms depending on the density and/or the number of atoms in the cold MOT cloud are well known 
' even for degenerate groundstate atoms 12J. They essentially rely on reabsorption of scattered photons in the cold 
^ MOT 13]. This situation is encountered as soon as light multiple scattering sets in, i.e. when the light scattering 
^ mean free path £ is comparable to the MOT size L. Similar multi-atoms and multi-photons effects certainly also exist 
^ with earth-alcaline MOTs, but they can be made negligible by working at low densities (n « 2.5 lO^atoms/cm'^) and 
at low optical thickness (6 = L/£ « 0.3 at resonance). 

Strictly speaking, Doppler theory can also be tested with degenerate atoms provided a suitable molasses polarization 
configuration is chosen 14]. Indeed, for ID-cooling with parallel polarizations ((7+ — (t+ or tt — tt configurations), no 
Sisyphus cooling can occur. However, even in this case, the authors found unusual high temperatures. They suggested 
that this discrepancy was due to molasses intensity imbalance, leading to a local drift of the atomic average velocity 
in the cloud. We agree with this explanation as we will show in detail in this paper. 

We have performed specific temperature measurements on ^^Sr (Sr in short hereafter) in a ID-cooling configuration. 
This is the ideal situation to test the Doppler theory for two reasons. First, we induce cooling on a = — > Je = 1 
atomic optical transition. Second, 2D or 3D Doppler cooling would be more difficult to analyze quantitatively. Indeed 
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interference between the six laser MOT beams can induce light shifts modulations, leading to modifications of the 
usual friction and diffusion coefficients |15|. These spurious effects do not exist in fD configuration with mutual 
orthogonal polarized counter-propagating laser beams. The new important ingredient here is that we incorporate in 
the original Doppler theory the effect of transverse spatial intensity fluctuations of the laser molasses profiles. 

The paper is organized as follows: after some details about the experimental procedure, we compare velocity 
dispersion measurements with the Doppler theory predictions f section Hljl . We have found temperatures 10 times 
larger than predicted. Indeed, at a laser detuning 6 — — r/2, we have measured a velocity dispersion ay = 0. 7±0. 2 m/s 
(T « 5mK) whereas the Doppler theory predicts — 0.23m/ s [T k, 0.5mK). In section ITlll we will show, on very 
general arguments, why the molasses transverse intensity fluctuations heat the cold cloud. We then derive an analytical 
model valid in two limiting cases: when the transverse distance Lj_ travelled by an atom during the longitudinal 
velocity damping time is much shorter or much longer than the transverse correlation length of the molasses 
intensity fluctuations. As expected, these two different limits give rise to different final velocity distributions. In 
subsection IIIIDI we will compare our results to a Monte-Carlo simulation. Most of our experimental data correspond 
to L± ^ ^s- In this case, the dynamics of the velocity distribution shows an unusual behavior. The velocity 
dispersion is first reduced by the friction force and then increases after a time scale related to (section llVfl . 
This specific behavior confirms the dominant role played by the transverse intensity fluctuations in ID-cooling with 
a Jg = — > Je = 1 transition. 

II. EXPERIMENTAL RESULTS 
A. Magneto-optical trap 

The cold strontium cloud is produced in a MOT. The Jg = ^ Jg = 1 dipole transition under consideration 
is the optical atomic line ^5*0 — ^Pi at A = 461 nm. The excited-state natural linewidth is r/27r = 32 MHz and 
the corresponding saturation intensity is Is = 42.5 mW/cm^. First an effusive Sr beam is extracted from a 500° C 
oven. Then a 27 cm long Zeeman slower reduces the Sr longitudinal velocity within the velocity capture range of 
the MOT 50 m/s). The Zeeman slower, MOT, and probe laser beams operate at 461 nm and are generated 
from the same frequency-doubled source detailed in [3. Briefly, a single-mode grating stabilized diode laser and a 
tapered amplifler are used in a master-slave conflguration to produce 500 mW of light at 922 nm. This infrared light 
is then frequency-doubled in a semi-monolithic standing-wave cavity with an intra-cavity KNbOa nonlinear crystal. 
The cavity is resonant for the infrared light while the second harmonic exits through a dichroic mirror providing 
150 mW^ of tunable singlc-modc light, which is then frequency locked on the 461nTO Sr line in a heat pipe. We use 
acousto-optic modulators for subsequent amplitude and frequency variations. The MOT is made of six independent 
trapping beams. Each beam is carrying an intensity of 5.2mW/cm^ and each beam waist is 8 mm. The trapping 
beams are red-detuned hy S = —T with respect to the atomic resonance. Two anti-Helmoltz coils generate a 70 G/cm 
magnetic field gradient to trap the atoms. The number of trapped atoms, as deduced from fluorescence measurements, 
is iV « 2.5 10^. The (Gaussian-shaped) cloud dispersion is roughly 0.6 mm. The velocity dispersion of atoms in the 
MOT is typically Av ~ Im/s. Technical details about temperature measurements are given in the following section. 

B. ID-cooling and time-of- flight measurements 

In addition to the previously discussed experimental set-up (MOT, Zeeman slower), two additional contra- 
propagating laser beams are used to perform a ID molasses (see figure^. These ID cooling beams are located 
in the horizontal plane. Some quarter-wave and half-wave plates are used to fix the relative polarizations of the two 
contra-propagating laser beams. An acousto-optical modulator in a double-pass conflguration is used to adjust the 
laser frequency from resonance down to — F. The same acousto-optical modulator is used to control the laser intensity 
up to Is. The beam waist is 3.5 mm at the cold cloud position, thus much bigger than the MOT size. 

The time sequence of the experiment is generated by a personal computer with digital output ports dedicated to 
this task. An internal clock updates and synchronizes the digital output ports every 10 /xs. This elementary time 
step is short enough for our purposes. The digital output ports are connected to switches which turn on and off the 
lasers beams (rising and falling times shorter than 1 /zs), the magnetic field (rising time of few ms and falling time of 
100 fis) and the CCD camera chopper (opening and closing times of 200 /is). The time sequence is designed as follows: 
first the MOT is operated during about 20 ms. Then the MOT lasers and magnetic gradient are switched off. The 
ID-molasses laser beams are then switched on during 500 /is. This cooling time is appropriately chosen : atoms reach 
the cooling steady-state regime while, at the same time, the expanding atomic cloud remains smaller than the cooling 
beams size. After the ID-cooling sequence, the temperature of the cold cloud is extracted from a time-of-flight (TOF) 
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measurement technique. For this purpose, all the laser beams are switched off and the cold cloud expands ballistically 
in the dark. The duration of the dark period is varied from a few 100 fis up to 1.5 ms. Then a fluorescence image of 
the expanding cloud is recorded on a CCD camera by briefly switching on the MOT beams during 20 fis. The whole 
time sequence is then repeated as long as necessary to obtain a good signal-to- noise ratio (see figure^. 

Figure 13 shows three images of the cloud after different ballistic expansion times. The elliptical shape of the cloud 
is a signature of the ID-cooling sequence. Indeed the cooling axis is precisely the small axis of the ellipse. This 
dimension will be thereafter referenced as the longitudinal axis. In the two other dimensions, thereafter referenced 
as the transverse axes, the cloud is heated by random spontaneous emission. From these images, we extract the 
dispersion of the longitudinal spatial distribution of atoms in the cloud, integrated over the transverse directions. 
Because of the finite size of the image, we cannot have access to the long tail of the spatial distribution. Hence, 
we truncate all recorded distributions at 2% of their maximum value. Assuming the initial atomic positions and 
velocities in the MOT to be uncorrelated and centered, the ballistic time evolution of the spatial dispersion is simply 
cr^(i) — cr^(O) + a^t"^. The velocity dispersion = a/ (w^) is then easily extracted from the experimental data. To 
cross-check these measurements, we have also used a different measurement method, namely a spectroscopic technique. 
It consists in probing the cold cloud with an ultra-stable laser beam tuned on resonance with the ^Sq—^Pi Sr transition 
line (note that this is also a — s- 1 transition). Because this line is spin- forbidden, its frequency width is very small 
(7.5 kHz) and is thus Doppler-broadened in the MOT Hence, we have a direct access to the velocity distribution 
by measuring the Doppler-induced spectral width. We have checked that these two different techniques give the same 
results, thereby confirming the validity of the TOF measurements, easier to handle, and for this reason routinely used 
in this experiment for temperature measurements. 

C. Doppler theory and results 

Doppler theory is based on absorption-fluorescence cycles which induce a cooling mechanism competing with a 
heating mechanism. Cooling is due to a mean friction force F = —m'jv which damps the velocity and heating is 
due to a Langevin fluctuating force F^, (photon noise) giving rise to diffusion with constant D in velocity space. At 
equilibrium these processes exactly balance (fluctuation-dissipation theorem) and the Doppler temperature is found 
to be ksTu — ma^ — mD/j. All experimental data presented in this paper have been obtained in the lin\\lin 
polarization configuration. In this case, analytical expressions for the damping coefficient 7 can be found in and 
in for the diffusion constant D. The analytical expression of the Doppler variance at low laser intensities / and 
negative detuning S is then 

Here A = 26 /T is the laser detuning in units of r/2 and sq = I /Is is the on- resonance saturation parameter. As 
could have been guessed, the Doppler temperature is related to the only energy scale of the problem, namely the 
excited-state energy width TiT. We thus see that, at low intensities, ao ~ o'o(A) {1 + (3{A)so) grows linearly in sq with 
slope /3(A)(To(A) from (To(A). Minimization of o'o(A) is achieved for A = — 1 and gives /3 — 3/7, (Tq — y^ThV /20m « 
0.23 m/s. 

In the strict ID Doppler theory, the final velocity distribution and the corresponding temperature should depend 
on the polarization configuration. For example, in the lin\\lin polarization channel, photon redistribution processes 
between the two contra-propagating laser fields can occur, whereas they are forbidden in the lin _L lin polarization 
channel. Hence, the polarization configuration affects both friction and diffusion. We have tested different polarization 
configurations {lin\\lin, lin _L lin and cr"*" — a~) but no significative temperature modifications have been found in the 
parameters range used in this experiment (see below). 

FigureElshows the experimental variation of the velocity dispersion as a function of A at low intensity (so = 0.08). 
The global behavior is the same as the one predicted by the Doppler theory but systematically higher. For example, 
at A = —1, we get cr„ — 0.7 ± 0.2 m/s {T « 5mK) whereas the Doppler theory predicts di, w o-q ~ 0.23 m/s 
(T sa 0.5 mK). The mismatch is even more pronounced at small detuning and tends to be reduced at high detuning. 
Figure^] obtained at A = —1, shows a strong linear increase of ay as a function of sq. The measured slope is 0.9 m/s 
whereas the Doppler theory only predicts 0.1 m/s. 

The on-resonance optical thickness corresponding to the data in figures 13 and 01 is typically b « 0.3 when the 
experimental sequence starts. At the end of the sequence, because the atomic cloud expands, b ~ 0.15. These optical 
thicknesses are not sufficiently low to discard multiple scattering effects (this would require 6 ^ 1). The role of 
multiple scattering of light is two-fold. First, it induces photon reabsorption leading to an average repulsion force 
and heating. However at 6 w 0.3, heating, as observed in can here be neglected. Second, it implies beam 
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attenuation (Lambert-Beer law). This means that the average optical force is weaker for atoms located deep in the 
cloud and we get, as a net effect, an average compression force. In a MOT, these compression and repulsion forces are 
equilibrated by the trapping force This is not the case in optical molasses. In a ID-molasses, the compression 
force is expected to dominate over the repulsion force induced by photon reabsorption, at least at moderate optical 
thicknesses. A quantitative estimate of multiple scattering effects is thus not easy. For example, in our case, the 
Lambert-Beer law predicts a maximum relative intensity imbalance of approximately 23% at low intensity, which is 
not negligible. However our experimental results strongly suggests that multiple scattering velocity inhomogeneous 
broadening is small. Indeed, as sq is increased, the optical thickness decreases since the scattering cross-section is 
reduced. The compression force is thus decreased and cannot explain the strong increase of vs sq evidenced in 
figure^ Furthermore, multiple scattering induces some correlations between position and velocity which should alter 
the cloud ballistic expansion. This has not been observed in figurc[21 As an ultimate test, we have change the number 
of atoms in the MOT by a factor 3 without detecting any modification of a^. 

III. COOLING WITH INTENSITY IMBALANCE 

We analyze in this section how spatial stationary transverse intensity fluctuations can modify the Doppler cooling 
predictions in a quantitative way. Starting from the analytical expression of the average force in ID molasses, we 
take into account these spatial intensity fluctuations and we derive analytical results in two limiting cases, namely 
when cooling is achieved before intensity fluctuates and in the opposite case. This analytical model neglects the 
photon noise encapsulated in the Langevin force F^, leading to a nonzero Doppler temperature. However, for a 
quantitative comparison with experimental data, we have developed a Monte-Carlo simulation which fully takes into 
account all fluctuating mechanisms (molasses intensity fluctuations and photon noise). These results will be detailed 
in subsections IIII Cl and IIII Dl 

A. Origin of spatial intensity imbalance 

In the standard Doppler theory, the molasses beams are described as perfect plane waves. In real experiments 
however, the beams have a Gaussian-shaped transverse profile. Usually one can ignore the transverse dimensions of 
the beams because they are generally much larger than the MOT size. Nevertheless, the transverse intensity profiles 
are not defect-free ideal Gaussian profiles. Indeed, even starting from diffraction-limited laser beams, imperfections 
of optics elements (dust, aberrations, etc) induce scattering which generates an intensity speckle. Of course, in well- 
controlled experiments, this speckle pattern remains relatively small compared to the average beam intensity, but 
those two fields add coherently. In our experiment, we have measured the laser beam spatial transverse fluctuations 
by placing a CCD camera at the approximative position of the MOT. Subtracting the ideal Gaussian profile, we have 
computed the intensity fluctuations histogram (see figurelSJ. We found a Gaussian histogram with a relative standard 
deviation (with respect to the average intensity) in the range 10 — 20%. 

The ID-molasses is created by two contra-propagating beams issued from the same source. However, the speckle 
intensity generated by the optics imperfections in each arm are independent. It is thus reasonable to consider, at 
each transverse position in the beam profile, the intensities li {i = 1,2) of the two molasses beams to be random 
independent Gaussian variables. Noting by Si the corresponding on-resonance saturation parameters, we assume the 
probability distribution V{si) to be the same for the two beams. Hence the (common) first two moments, at each 
transverse position, are (s^) — sq and cr^ = ((si — sq)^) (i = 1,2). In the following wc will characterize the molasses 
intensity fluctuations by the ratio Ts = cFs/'^sq. From experimental data, we have in the range 5 — 10%. Another 
important feature is the transverse spatial correlation function (assumed to be the same for each beam), namely 
C(r) = (si(r')si(r' -t- r)) (i = 1,2). The characteristic decay length of C(r) defines the transverse correlation length 
^s- Strictly speaking, because of diffraction effects, cannot be zero. We have experimentally found to be of the 
order of few tenth of (^s = 30/xm for the example shows in figure (SJ. 

B. Analytical model 

For a ^ 1 transition in the lin\\lin configuration, atoms are modelled by a two-level system involving the same 
Zeeman states for each molasses beam. The two beams have the same red-detuned frequency (J < 0), opposite wave- 
vector k but locally different saturations which fluctuate independently across the transverse profile of the beams 
but with the same correlation length • The mean local radiation pressure force experienced by a single atom, with 
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longitudinal velocity v, takes the simple following form: 
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where vr = fik/m is the recoil velocity Qmrn/s for Sr). We have here assumed that the mechanical actions of 
the two laser beams can be added independently, only considering the total saturation of the optical transition to be 
(si + 82)- We have furthermore replaced (si + S2) by 2so in the denominator. A more careful analysis, along the lines 
given below, shows that this is indeed valid with our experimental parameters. 

Near the steady state, the Doppler broadening become negligible, kv ^ (F, 5) and the force can be safely approxi- 
mated by: 



F « — to(si + S2)7„ w + to(si — S2) a 
w —2msQ^yV + m{si — S2)a 

Defining the recoil angular frequency ljjr, — 2ttvr = hk'^/2ni {vr, k, IQ.QkHz for Sr), we have 
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The friction force is given by the term proportional to v and gives rise to a mean damping time: 

Ty = {2S0"fv)~^ 



(4) 
(5) 



(6) 



(7) 



Consider now an atom moving across the molasses beams with transverse velocity vj_ . This velocity is unaffected by 
the longitudinal cooling. So, as the atom flies across the beam, it may experience induced temporal molasses intensity 
fluctuations. The induced correlation time of these fluctuations is simply 



^s/vjL 



(8) 



We have thus two competing dynamical processes, the cooling one with characteristic time t„ and the intensity random 
(induced) temporal variations with correlation time r^. The global velocity dynamics can be easily analyzed when 
these time scales are well separated. 



1. The regime r„ ^ Ts 



For atoms fulfilling this condition, the damping process is completed before intensities fluctuate. Then the constant 
force and the friction coefficient in Q can be considered as time-independent. Atoms thus behave as if v± was 
zero. However, because is much shorter than the cloud size, each atom experiences randomly fixed distributed 
intensity imbalances and, in turn, a random constant force. The random stationary velocity is Voa = aTy{si — S2). It 
is Gaussian-distributed with zero mean and dispersion 

/- /- F 1 + |A|2 + 2so 

(Too = v2ar„ (Ts = V2 — (9) 

k 4| A| 

To find the final velocity distribution, one has to incorporate the photon noise force F^. As the photon noise and 
intensity fluctuations are independent, the total velocity variance is the sum of the Doppler and intensity variances 
'^v — ~^ '^oo- this regime, intensity fluctuations thus give rise to an extra-heating mechanism inducing an 
inhomogeneous broadening of the Doppler velocity distribution. At A — —1, sq small and r^ — 9%, one gets 
aoo ~ 0.94 m/s and ~ 0.97 m/s. 



2. The regime Tv 3> Ta 



In this regime, the molasses intensities fluctuate wildly before the atom reaches its stationary state. Thus they 
cannot be treated as time- independent quantities. They can however be treated as independent Markovian (short 
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memory) processes. From the expression (0} of the total force, we derive the following master equation for the velocity 
distribution ?<(«) (details will be given elsewhere): 

P{v,t + Ts)^ J dsi dS2 dv'V{si) V{S2) P(w', t)5(^~v' ~ ^T^j (10) 

In this equation, at each time step Ts, the saturations si and S2 take on new values uncorrelated with the previous 
ones. 

One should note that we have kept here both the additive noise, given by the (si — S2) term, and the multiplicative 
noise, given by the (si + S2)v term. This last term had been replaced by 2sqv in our previous analysis (expression 
(O of the force). In the limit of small velocity changes at the elementary time scale r^, one can derive from (|10|l a 
Fokker-Planck type equation (see Appendix IVII|) : 

^ = ^|;(^P(M)) + |;(i^(.)P(M)) (11) 
Defining D^c = o'^/2ti,, the velocity-dependent diffusion constant has the following expression: 

D{v) =D^[\ + {1 + 2r2) [v/aj" ]Zl^D^[l + [v/a^f] ^ (12) 

since generally 2r^ <C 1. This velocity-dependence originates from the multiplicative noise and leads to abnormal 
diffusion. The stationary velocity distribution Pq is then easily found to take the following form (M is a normalization 
constant): 

¥o{v)^U[l + {v/cr^f]-(^+^) (13) 
where rj = Tv/ts ^ 1. The full width at half maximum of this distribution is 

< (Too (14) 

As easily checked D[v) ~ Tg/r^ for v ^ T^. Hence the impact of abnormal diffusion is in fact negligible when 
Tv ^ Ts and the velocity distribution is very close to a Gaussian with velocity dispersion 

< o-oo (15) 



which is much narrower than the velocity dispersion found in the long correlation limit (^J by the factor yTs/ry. 

Figure El shows as a function of in the long and short correlation time limits (A = —1, sq ^ !)■ When ^ 0, 
(T„ goes to zero in both limits. This does not mean that the temperature goes to zero in experiments. Indeed, in 
the previous discussion we have discarded the photon noise at the physical origin of the bare Doppler theory. Its 
contribution at low saturation is given by the dotted line in figure El We can see that the noise induced by long 
correlation times starts to dominate at > 2.5%, a fairly small value. Hence, to achieve cooling up to the Doppler 
limit at A = —1, one has to minimize the ratio Tc/t„ for a fixed value of r^. This can be done in different ways. First, 
by decreasing the laser intensity to increase t^. Second, by decreasing Ts. As shown by ((HJ, this can be achieved 
either by avoiding large transverse intensity defects in the laser profiles to decrease the correlation length or by 
increasing v±. This is indeed was has been observed in atomic beam experiment |2l| where temperatures very close 
to the Doppler limit had been found along one transverse dimension. 

In a MOT, if the fluctuations of the transverse intensities are not the same for all three dimensions, the temperature is 
expected to be anisotropic. Moreover coupling mechanisms occur between the cooling dimensions. Thus high velocity 
dispersion along one dimension tends to reduce the temperature in the orthogonal plane. In our ID-molasses, the 
transverse velocity is fixed by the initial MOT sequence and Ts/t^ is usually bigger than 1. As we will see in Section 
III.D., this means that transverse intensity fluctuations are the major heating mechanism. 



C. Monte-Carlo simulations 



In order to quantitatively test our theory and not relying on questionable approximations, we have developed a 
Monte-Carlo (MC) simulation. This MC simulation fully takes into account the photon noise leading to the Doppler 
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cooling limit, saturation of the transitions and the transverse intensity fluctuations discussed so far, i.e. with arbitrary 
correlation time Tg and saturation dispersion (Tg. 

The shortest time scale in the MC calculation is given by the excited-state lifetime Tg — l/T. To mimic the 
transverse flight of atoms, the molasses saturation fluctuations Ssi = {si — sq) (« — 1,2) become time-dependent 
parameters and evolve according to the damped random path discrete equation (n labels the number of time steps 
re) 

6s,{n+l) = (l-p)Ss,{n) + Rn (16) 

with p <C 1. The first term in this equation is the friction term, relaxing saturation to its stationary value so, while 
the last term is a random variable, with zero mean value, uniformly distributed over the range [—e/2, e/2]. A simple 
calculation shows that (i?^) = e^/12. We further assume that the i?„'s are decorrelated. The continuous limit of Hlt)|) 
is 

^^ + -^Ss,^R{t) (17) 

at Te 

As the correlation time should correspond to the damping time of this equation, we see that p = T^/Tg. Thus fixing 
Tg fixes p in l|l()|l . The last term of (|17|) is a (5-correlated Langevin term with zero mean. Its time-correlation function 
is {R{t')R{t)) ~ 2'DS{t — t') where it is easily shown that V = e^/24Te. The fluctuation-dissipation theorem then 
dictates al = VTg, leading to 



t^2y/&pag (18) 

This last result can be found more elegantly by squaring (|16|l and averaging over the probability distribution of i?„. 
Hence experimental determination (or convenient choice) of the macroscopic ingredients Tg and ag fixes in principle 
the microscopic ingredients e and p in (|15|l . 

Figure shows the obtained final velocity dispersion as a function of Tg/r^ for sq = 0.04 and = 9%. The 
velocity damping time, calculated with 10, is Tt, = 200 /xs. As expected, (Tv is higher for long correlation times 
Tg ^ Ty. When Tg ^ Ty, ay is very sensitive to Tg. At lower values Tg <JC t„, <7y is minimum and, for these parameters, 
reaches the Doppler limit. The final distributions (not shown here) are quasi-Gaussian, even at small Tg. 



D. Quantitative comparison with experimental data 



We have showed in section II. C that the measured velocity dispersion were always larger than the Doppler theory 
predictions. We now compare our experimental data with the results of our previous MC simulation. In order to 
stick as close as possible to the experiments, we also take into account the transverse velocity distribution (a centered 
Gaussian with dispersion Av± « 0.8m/ s) in the MC simulation. 

In figure O the solid curve corresponds to the MC simulation performed at = 9% and = 60 pm. These 
quantities are fit-parameters in the MC simulation but remain in the range of the measured ones fsee llll A^ . As one 
can see, the agreement with experimental points is now very good, providing a clear understanding of the physics 
at work in the experiment. Coming back to lO, close to resonance, the friction term decreases while the intensity 
imbalance term increases. Hence, the mismatch between experiment and Doppler theory is maximal. 

Again, an excellent agreement between MC simulation (solid curve) and data is found in figure ^ To properly 
understand these results, one has to remind that the cross-over region between the long- and short-correlation time 
limits is given by Tg — Ty, or equivalently by sq — v±/2^y£^g ~ 0.1. TOF measurements before the molasses sequence 
have shown that v± « 0.8 m/s whereas £,g is the chosen MC parameter. This is exactly what is seen in figure0]and the 
Doppler theory is recovered when sq 0. However, this short-correlation time limit (r^ <C Ty) was not experimentally 
accessible because the ID-cooling duration sequence was not long enough to reach the steady state. This is why most 
of the data points correspond to the long-correlation time limit (r^ ^ Ty) and can be compared to ©. This is done in 
figure 21 where the dashed line corresponds to the prediction given by ©. The general behavior is correct, indicating 
that saturation of the atomic transition indeed plays a significant role. However the prediction is a little bit too 
large. This is not surprising since, for the explored range of parameters, some atoms will always have sufficiently high 
transverse velocities to fulfill the short correlation time criterion, both in the experiment and in the MC simulation. 
As a consequence, the velocity dispersion ay will be reduced. In other words, © corresponds to a zero-transverse 
velocity case giving over-estimated predictions. 

Figure IHl shows the longitudinal spatial distribution obtained for the longest ballistic time at A = — 1 and sq = 0.08. 
It is two times broader than the initial spatial distribution and essentially proportional to the velocity distribution. 
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We can then try to compare it to the velocity distribution obtained with the MC simulation (solid line). We have also 
plotted the Gaussian distribution expected from standard Doppler theory (dashed line). The actual non-Gaussian 
shape of the MC distribution is explained by the transverse velocity dispersion (At;j^ w O.Sm/s). For each fixed 
transverse velocity, the distribution is quasi- Gaussian with a width depending on v±_ (see figure and discussions 
in Section IlIIBjl . For this simple reason, summing over the transverse velocity distribution leads to a non-Gaussian 
distribution. The non-Gaussian shape of the MC simulation matches the experimental distribution better than the 
Gaussian one. 

As a conclusion of this section, we again stress that we have quantitatively explained both the behavior of the 
velocity dispersion as a function of laser detuning and intensity and the observed non-Gaussian distributions. This 
puts strong evidence on the fundamental role played by molasses intensity defects in the cooling process. 

IV. COOLING DYNAMICS 

Figure El show the time evolution of cr„ for three different values of corresponding to the short, intermediate 
and long correlation time limit. When Ts is short (or equivalently v± large), the corresponding curve (a) displays an 
exponential-type behavior decaying to the Doppler steady-state value. This curve is in agreement with the dashed line 
Doppler prediction. When is long (or equivalently v± small), the dynamics evidenced by curve (c) is more complex. 
First ffy reaches a minimum value within a time scale corresponding to the damping time Ty . Then, cr„ increases again 
and reaches the steady state value Uoo predicted by @. When the mean velocity and the local intensity imbalance 
arc uncorrelated in the initial state, then this long correlation time behavior is generic and does not depend anymore 
on the initial state. Because the intensity correlation length is large, the heating mechanism takes also some time 
to build up and the initial velocity distribution starts first to shrink. This is easily explained by considering the 
dynamics induced by Q for times shorter than where transverse intensity fluctuations are dynamically frozen. 
After a proper averaging over intensities fluctuations, the following analytic expression for ay{t) is derived: 

a„(t)2 = (Aw2+a^)exp(-2VT.)+a^(l-exp(-Vr„)) (19) 

where Aw corresponds to the initial longitudinal velocity dispersion, the small time expansion t ^ Tt, of this equation 
gives: 

ay{t) ^ Av{l- t/Ty) (20) 

clearly evidencing the velocity spread narrowing at small times. We have experimentally tested this specific behavior. 
Once the MOT is loaded, we have switched off the cooling laser beams within a time window of variable duration 
Tdark- This dark sequence must be long enough to allow atoms to travel over transverse distances larger than the 
correlation length of intensity fluctuations. Then, once the cooling lasers are switched on again, any correlation 
between the mean atomic velocity and the local intensity imbalance is wiped out. We should thus observe the behavior 
predicted by curve (c) in figure |^ The velocity dispersion is measured after a time r by the TOF technique (see 
section^. In figure ITUT a). we plot ay as a function of r after a dark period Tdark = 0.5 ms. We reproduce nicely the 
corresponding theoretical prediction. First we observe a decrease of ay followed by an increase, up to the final value 
(which is here the same as the initial value). We also checked that the dark period has to be long enough, as evidenced 
in figure ITUT b). When Tdark is short enough (less than Ims), decorrelation between the mean atomic velocity and 
local intensity imbalances is not completed. When Tdark is long enough (larger than 1 ms), complete decorrelation is 
achieved and the dark period no longer plays any role in the cooling dynamics. 

For practical reasons, these experiments were done with the MOT {i.e. on an 3D-cooling configuration in the 
presence of the magnetic field gradient) whereas the MC simulations were done for a ID-cooling scheme. This major 
difference does not allow for a quantitative comparison between theory and experiment. However, the qualitative 
agreement is pretty good. This result suggests that the dominant extra-heating mechanism in a MOT is the same as 
in ID-molasses. 

Because the velocity distribution always starts by a compression period, one could imagine a cooling strategy 
implementing suitable repetitions of dark time windows to achieve Doppler-limited cooling. This simple idea however 
is not easy to handle because complete decorrelation between initial atomic velocities and intensity imbalances during 
these dark windows is requested. As the cloud cools down, the duration of the subsequent dark periods has to be 
increased accordingly to maintain this decorrelation. Unfortunately, as soon as the cooling time becomes very long, 
spurious effects such as large cloud expansions, then have time to set in. 
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V. CONCLUSIONS 

We have evidenced in this paper the important role played by transverse spatial intensity fluctuations in ID laser 
cooling of zero spin groundstate atoms where no Sisyphus cooling can occur. For intensity imbalanced molasses beams, 
the total radiation pressure force decomposes, at small velocities, into a friction force and a constant force. The latter 
is at the root of an additional heating mechanism. These two forces essentially depend linearly on the molasses 
intensities. In the presence of sub-Doppler cooling mechanisms, the friction term becomes intensity-independent 
whereas the constant force remains proportional to the intensity 0,|23il23- Hence, with sub-Doppler cooling, the 
effect of intensity imbalances can be arbitrarily small and, in turn, does not play any significant role. This is completely 
different for Doppler cooling where this effect remains dominant in most experimental cases. 

The cooling steady-state reached by the atoms depends sensitively on the ratio C, between the correlation length of 
transverse intensity fluctuations and the transverse distance travelled by the atoms before reaching the steady-state. 
For small C, the transverse fluctuations can be modelled by a Langevin force inducing an additional heating mechanism. 
For large C, the equilibrium state is reached at frozen molasses intensities. Atoms at different transverse positions 
then probe all possible intensities imbalance. This averaging procedure mainly affects the final atomic velocity and 
induces an inhomogeneous broadening of the Doppler velocity distribution. The impact on the final temperature is 
more severe than in the small C regime. This is evidenced by our experimental data which mostly lie in the large ^ 
regime. It is however possible, in principle, to reach the small C regime by appropriately reducing the intensities of 
the cooling beams and get final temperatures closer to the Doppler limit. 

The cooling dynamics also exhibits an unusual behavior in the large ^ regime : the atoms are first cooled down 
before being heated up. This specific feature is also found in our experiments, thus confirming the central role of 
transverse intensity fiuctuations in Doppler cooling. 
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VII. APPENDIX : DERIVATION OF THE MASTER EQUATION 

The Fokker- Planck type equation is obtained by starting with the master equation H1U|) : 

P(«, t + Ts) = I dsi ds2 d^;' V{si) V{S2) V{v', t) 6{v - v' - :^lfiifEl^^^) (21) 



By Fourier transforming this equation with respect to w, we get: 

P(q, t + Ts) = -^ I dsi ds2 Av' V{si)V{s2) PK, t) exp[-z<7(«' + ^^^^ll^^r,)] (22) 
V 27r J m 

We plug now expression I^J of the force in this equation to find: 

Piq,t + Ts) = [ dv"P[qT,ia-jW)] V*[qTs{a + -f.v')] Piv' ,t) exp{-iqv') (23) 



/2tt 

where the star denotes complex conjugation and where 



'P(-u) = / ds-p(s)exp(-ms) (24) 



is proportional to the Fourier transform of the molasses saturation distribution V. As discussed in section UlI Al the 
distribution 7^ is a Gaussian with mean sq and dispersion as- A Taylor expansion of (|23l) up to second order in 
then leads to: 



dP{q,t) 1 



/ dt;' Uqv' -q'' — [v'^{l + 2r'^) + 2r^^^] P{v',t) exp(-igw') (25) 
J L 2r„ 7^^ 



dt 

An inverse Fourier transform then gives equation once we note that ctqc = V^fsu/jy according to equation 
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FIG. 1: Schematic drawing of cnA nf +>io +;tt,o aom^onno ncorl in fVio ovnorimonf T^Vio oiv iri dependent MOT laser 

beams are along the x axis and isional cooling beams 

contrapropagate along the y ax integrated over z. 




[Time (ms)] 

FIG. 2: After the ID-cooling sequence along axis y, the spatial distribution in the horizontal plane is collected at three different 
ballistic expansion times by the time-of-flight technique. From these images, we extract the velocity dispersion av in the cloud 
(see text). 
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FIG. 4: Velocity dispersion Uv as a function of the mean saturation parameter so = I /Is at A = 25 /T — —1. The experimental 
data (circles) are compared to the Doppler prediction (dotted line) and to the Monte-Carlo simulation (solid line) at Vs = 9% 
and = 60 /xm. Whereas the Doppler theory is completely off, very good agreement is found with our theoretical model. The 
dashed line prediction corresponds to equation © (see text). 
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FIG. 5: Figure (a) shows typical relative intensity fluctuations of a laser beam once the smooth Gaussian profile is removed. 
The typical intensity dispersion is roughly 10% of the total signal. Figure (b) shows the spatial fluctuations correlation function 
C{x) vs a transverse coordinate. The transverse distance at which this correlation function vanishes defines the correlation 
length ^s. Experimentally we found — SO nm (diamonds). The oscillation on top of the decay is due to the interference 
generated by the front window of the CCD camera. The fast decrease observed at short distances is also present when the laser 
is off (crosses). This fast deer — — ;„ +u„ r<r<-r, — ™ 1 — -""ponds to the pixel size. 
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FIG. 6: velocity dispersion av for the long (solid line) and short (dashed line) correlation time limit as a function of the relative 
fluctuation of the intensity Ts = as/2so- In the short correlation time limit t, = SOfis. The dotted line corresponds to the 
Doppler limit at low saturation and A = — 1 (see text). 
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FIG. 7: Final velocity dispersion cr„ as a function of the correlation time Ts (full circles). The transverse and longitudinal velocity 
dispersions are both equal to Aw = 0.8 m/s. The mean saturation of each beam is so = 0.04 and the saturation fluctuation 
parameter is Ts ~ 9%. The upper dashed line represents the long correlation time limit while the lower one represents the 
Doppler-cooling limit (see textl. 




Velocity (m/s) 

FIG. 8: Velocity distribution obtained at so = 0.08 and A = —1. Circles : experiment. Solid line : Gaussian fit. Dashed line : 
Monte-Carlo simulation with rs = 7.5% and Ta = 20 fxs (see text). 




Time (|j,s) 

FIG. 9: Monte-Carlo simulations of the time evolution of the velocity dispersion ay for three characteristic correlation times: (a) 
short correlation time (ts = 1.25 /us) ; (b) intermediate correlation time (ts = 125 /xs) ; (c) long correlation time (ts = 1250 /js). 
The transverse velocity dispersion is Av± — 0.8m/ s, the mean saturation per beam is sq = 0.04 and the saturation fluctuation 
parameter is rs = 9%. The dashed exponential decay is the bare Doppier prediction (see text). 
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FIG. 10: (a) : Time evolution of the velocity dispersion av after a dark time window of duration Tdark ~ 0.5 ms ; (b) MOT 
velocity dispersion ay as a function of the dark window duration. The dark window starts 0.5 ms after the beginning of the 
cooling sequence (see text). 



